Private Sub ClearOutValues_Click () 

AtPOut.Text = an 
OutDeosity.Text = 
PwtOut.Text = no 
OutWehDepth . Text = » » 
OutN = 

OutBuniRate.Text = "" 
End Sub 



Private Sub ConvertButtan_Click () 

Dim BSModifier As Double 
Dim ti As Single 

Dim praax As Single, fmax As Single 
Dim pti As Double 

' Set blank input or output variables 
If FillBlankValuesO = False Then 

Exit sub 
End If 

' check if -input data loaded 
-If InLoaded = False Then 

MsgBox "Your need to load a data file." 
Exit Sub 
End If 

Turn Hourglass on 
MainModel .MousePoint er = vbHourglass 

Calculate throat erosion profile to match initial and final throat values 
If Thrust FromAt = o Then 
Call SetErosion 
-End If 

' Adjust reference burn rate to match theoretical web burned 
Call SetBuxnRate 

' Adjust C* until propellant weight burned matches given weight 
Call SetCStar 

' Plot input bum back profile 

Da taCutFrame. Visible = False 
StatusLabel. visible & True 

StatusText .Caption = "Plotting Humback Profile..." 
StatusText .Visible = True 




ProgressBarl. Value = 0 
ProgressBarl. Visible = True 
WehAtMaxPr ess. Caption = " n 
Refresh 



1.4 



MainModel .ThrustGraph .ScaleHeight = MaxThrust 
MainModel .ThrustGraph. ScaleWidth = MaxTime * 1.4 
MainModel .ThrustGraph. ScaleTop = MaxThrust * 0.1 
MainModel .ThrustGraph. seal eLe ft « -MaxTime * o.l 
MainModel .BBPGraph. ScaleHeight = MaxBS * 1.4 
MainModel . BBPGraph. ScaleWidth = MaxWeb * 1.4 
MainModel .BBPGraph. Seal eTop = MaxBS * 0.1 
MainModel .BBPGraph. ScaleLeft = -MaxWeb * 0.1 

BBPGraph.Line (0, MaxBS * 1.3) - (0, MaxBS * 0.1) 
BBPGraph.Line (0, MaxBS * 1.3) -(MaxWeb * 1.3, MaxBS * l 
-For i = 2 To InPoints 

r rf i / Int (InPoints / 10) = Int(i / Int (InPoints / 10)) Then 
ProgressBarl. Value = i / InPoints * 100 
Refresh 
End If 

BBPGraph.Line (WBurnln(i - l) , MaxBS * 1.3 - BSurflnfi 
* 1.3 - BSurfln(i)}, RED 
-Next i 

BBPGraph.Line (WBurnlnfi - 1), MaxBS * 1.3 - BSurflnU - 
MaxBS * 1.3), RED 



•3) 



1) ) - (WBurnln(i) , MaxBS 



1) ) - (WBurnIn(i - 1) , 



StatusLabel. Visible = False 
StatusText .Visible = False 
ProgressBarl .visible = False 
DataCut Frame. visible = True 
Refresh 



1 Write -burn back file to disk 
If BBPCheck. Value « checked Then 
fnum = FreeFile 

BbpFileSpec « BBPFileName .Text 
On Error Resume Next 
Open BbpFileSpec For Output As #fnum 
If Err Then 
MsgBox Error $ 
Exit Sub 
L End If 
Print #fnum, 
Print #fnum, 
Print #fnum. 
Print #fnum, 
Print #fnum, 
Print #fnum f 
Print #fnum, 
Print #fnum. 
Print #fnum f 

Chr$(9) & "Cf" t Chr$(9) 
Print #fnum, "" 



"Source: 11 & Chr$(9) 
"Temp: " & Chr$ (9) & 



& InFileSpec 
InTemp.Text 
Chr$(9) & AtO.Text 



"Initial Throat:" & 

"Final Throat:" & Chr$(9) & AtF.Text 
"Prop Wt:" & Chr$(9) & Pwt.Text 

"Burn Rate ® lSOOpsi:" & Chr$(9) & inBurnRate.Text 
"C*:" & Chr$(9) & CStar.Text 

n ii 

"Web Burned" & Chr${9) t "Burn Surface" & Chr$(9) & "Throat Area" * 



For i = 1 To InPoints 

Print #fnum, WBurnin(i) & chr$(9) & BSurfln(i) & chr$(9) & Atln(i) & chr$(9) 
& cfln(i) 
•-Next i 



Close fnum 
End If 



FIG. 2B 



' Generate converted file 
-If Progr adoptions (l) = True Then 
BSModifier « GenerateOutput ( } 



Plot Output thrust trace & burn back profile & write output file 



DataCutFrame. visible = False 

StatusLabel. Visible = True 

StatusText .Caption = "Plotting Output data. 

StatusText .Visible = True 

ProgressBarl .Value = 0 

ProgressBarl .Visible = True 

Refresh 



If OutputCheck. Value = (Checked Then 
fnum = FreeFile 

OutFileSpec = Out put FileName .Text 
On Error Resume Next 
Open OutFileSpec For Output As #fnum 
r lf Err Then 

MsgBox Error$, , "Error Opening Output File" 
Exit Sub 
L End If 
Print #fnum, 
Print #fnum, 
Print #fnum, 
Print #fnum, 
Print #fnum, 
Print #fnum. 
Print #fnum, 
Print #fnum, 
Text 

Print #fnum, 
Print #fnum, 
Print #fnum, 
Print #fnum, 



"Source:" & Chr$(9) & inFilespec 
"From Temp : 11 & chr${9) & inTemp.Text 
"To Tempi" & chr$(9) & Out Temp .Text 
"initial Throat:" & Chr$(9) & AtO.Text 
"Final Throat:" & Chr$(9) & AtF.Text 
"Prop wt:" Sl chr$(9) &. Pwt.Text 

"Burn Rate @ lSOOpsi & Input Temp:" & Chr$(9) & InBurnRate . Text 
"Burn Rate @ lSOOpsi & Output Tempi" & Chr$(9) & outBurnRate. 

"C*: M & Chr$(9) & CStar.Text 
"BSModifier" & Chr${9) & BSModifier 

n ii 



"Time" & Chr$(9) & "Pressure" & Chr$(9) & "Thrust" & 
"Throat Area" & Chr$(9) & n cf» & Chr$(9) & "Web Burned" & chr${9) 
Surface" 
Print #fnum, "" 
-End If 



Chr$(9) & 
& "Burn 



ThrustGraph . Cls 
BBPQraph . Cls 



BBPGraph.Line (0, MaxBS * 1.3) - (0, MaxBS * 0.1) 
BBPGraph . Line (0, MaxBS * 1.3)-(MaxWeb *" 1.3, MaxBS * 1.3) 
ThrustGraph. Line (0, MaxThrust * l.3)-(0, MaxThrust * 0.1) 
ThrustGraph. Line (0, MaxThrust * 1 . 3) - (MaxTime * 1.3, MaxThrust * 1.3) 

-For i = 2 To InPoints 

If i / int (InPoints / 10) = int(i / Int (InPoints / 10)) Then 
ProgressBarl. Value = i / InPoints * 100 
Refresh 
Lsnd If 



ti = ti + (ThrustOut (i) + ThrustOut(i -1)) / 2 * (TimeOut(i) - TimeOut{i - , 
1) ) * 
pti m pti + (PressOut(i) + PressOut(i - l) ) / 2 * (TimeOut(i) - Timeout (i - -, 
1) ) Tr 
r If Pressout(i) > pmax Then 
pmax = PressOut(i) 

WebAtPmax = WBurnOut(i) 2C 



End If 

If ThrustOut (i) > finax Then 

fmax = ThrustOut (1) 
End If 

If bsmax < BSurfOut(i) Then 

bsmax = BSurfOut(i) 
End If 



ThrustGraph.Line (Timelnd - 1), MaxThrust * 1.3 - Thrustlnd - 1) ) - (Timeln (i 
), MaxThrust * 1.3 - Thrustlnd)), RED 

ThrustGraph.Line (Timeout (i - l) , MaxThrust * 1.3 - ThrustOut d - 1) ) - 
(Timeout d) , MaxThrust * 1.3 - ThrustOut ( i )) , BLUE 

BBPGraph.Line (WBurnlnd - 1), MaxBS * 1.3 - BSurflnd - 1) ) - (WBurnln (i) , 
MaxBS * 1.3 - BSurfIn(l)) / RED 

BBPGraph.Line (WBurnOut d - 1} , MaxBS * 1.3 - BSurfOut d - 1) ) - (WBurnOut (!) , 
MaxBS * 1.3 - BSurfOut d) ) # BLUE 

-If Output Check. Value = Checked Then 

Print #fnum, Timeout (i) & Chr$(9) & PressOut(i) & Chr$(9) & ThrustOut (i) & 
Chr$(9) & AtOut(i) & Chr$(9) & CfOut d) & Chr$ (9) & WBurnOut (i) & Chr$ (9) & 
BSurfOut (i) 
End If 



L-Next i 



plf Out putCheck. Value = Checked Then 

Close fnurn 
1 End If 



ThrustGraph.Line (Timelnd - 1), MaxThrust * 1.3 - Thrustlnd - 1) ) - (Timelnd - 
1 ), MaxThrust * 1.3 - 0) , RED 

ThrustGraph.Line (TitneOut (i - 1), MaxThrust * 1.3 - ThrustOut (i - 1) ) - (TimeOut (i 
- 1), MaxThrust * 1.3 - 0), BLUE 

BBPGraph.Line (WBurnlnd - 1) f MaxBS * 1.3 - BSurflnd - 1 ))- (WBurnln ( i - l) , 
MaxBS * 1.3 - 0 ) , RED 

BBPGraph.Line (WBurnOut {i - l) , MaxBS * 1.3 - BSurfOut (i - 1 ))- (WBurnOut ( i - l) , 
MaxBS * 1.3 - 0), BLUE 

End If 

BurnTimeOut .Text = Format (Timeout (i - 1), "0.000") 
TIOut.Text - Format (ti, "0" ) 
PTIOut.Text = Format (pti, "0") 
PmaxOut.Text « Format (pmax, "0") 
FmaxOut .Text = Format (fmax, "0") 

WebAtMaxPress . Caption = Format (WebAtPmax , "0.00") 

StatusLabel .Visible = False 
StatusText .Visible = False 
ProgressBarl .Visible = False 
DataCut Frame. Visible = True 
Refresh 

' Turn Hourglass off 
MainModel .MousePointer = vbDefault 

End Sub 



Private Function FillBlankValues () As Integer 

FillBlankValues = False 

FIG. 2D 



If InTemp.Text = on Or OutTemp.Text - »■ Then 

MsgBox "You must enter both a Firing Temperature and an Output temperature." 

OutTemp . SetFocus 

Exit Function 
End If 

-If AtO.Text = »» Then 

MsgBox "You must enter an initial throat diameter . " 
AtO. SetFocus 
Exit Function 
End If 

If AtF.Text ■ »» Then 

MsgBox "You must enter a final throat diameter." 
AtF . SetFocus 
Exit Function 
-End If 
If Pwt.Text o «» Then 

MsgBox "You must enter a propellant weight." 
Pwt . SetFocus 
Exit Function 
-End If 
If InDensity.Text = "" Then 

InDensity.Text = Format (0 . 065, ".0000") 'eq(lO) 
-End If 

-If InWebDepth.Text = '»» Then 

MsgBox "You must enter a web thickness . ■ 

inWebDepth . SetFocus 

Exit Function 
End If 
-If InN.Text = "" Then 

MsgBox "You must enter a pressure exponent." 

inN. SetFocus 

Exit Function 
-End If 
If SigmaP.Text = D n Then 

MsgBox "You must enter an temperature sensitivity coefficient." 

SigmaP. Set Focus 

Exit Function 
-End If 

InBurnRate.Text « Format (0 . 0006 * Val (InTemp.Text) + 0.52, ".000") 

If Pwt Out. Text = " " Then 
PwtOut.Text o Pwt.Text 
-End If 

-If OutDensity .Text = "" Then 

OutDensity .Text = InDensity.Text 'eq(lO) 
-End If 

-If OutWebDepth.Text = "" Then 

OutWebDepth.Text = InWebDepth.Text 
End If 
-If OutN.Text = »" Then 
OutN.Text = InN.Text 
-End If 

-If AtFOut.Text = 11 ■ Then 
AtFOut.Text = AtF.Text 
-End If 

FillBlankValues = True 

End Function 



Private Sub Fonn_Load() 



FIG. 2E 



Private Sub SetErosionO 

'This subroutine varies the coefficient a in the equation E-aP^b 

'until the final throat calculated matches the given throat. 

*E ia the incremental radial throat erosion and P is the incremental 

'chamber pressure. The coefficent b is fixed at .924. This routine build the 

throat area area that becomes 

'part of the burn back file 

Dim ThisValue As Double, Goal As Double 

Dim Step As Double, Direction As Integer, Accuracy As Double 

Dim a As Double, b as Double, c As Double, d as Double, Trial As integer 

DataCutFrame. visible = False 
StatusLabel .Visible = True 

StatusText .Caption = "Calculating Erosion Coefficients..." 

StatusText .Visible = True 

ProgressBarl .Value = 0 

ProgressBarl. Visible ■ True 

Refresh 

'Calculate erosion coefficients 

a = 0.00001477 
b = 0.924 

Goal = val ( At F. Text) 
Accuracy = o.oooooooi 

Direction = 0 'first time thru step will be halved 

Step = 0.00001 'before it is applied 

Trial =0 

-Do 

Trial = Trial + 1 
ThisValue = Val (AtO .Text) 
EarlyEnd = InPoints 
For i = l To InPoints 

If i / Int (InPoints / 10) = lnt(i / Int (InPoints / 10)) Then 

StatusText .caption = "Calculating Erosion Coefficients... Trial #" & Trial 
ProgressBarl. value = 1 / InPoints * 100 
Refresh 
End If 

ThisValue = ThisValue + (a * Pressing) A b) *• (Tiirteln(i) - TimeinU - 1)) * 
2 'eq(2) 

Atln(i) - ThisValue A 2 / 4 * Pi 4 eq(3) 



FIG. 3.A 



If i > InPoints / 2 And Pressln(i) < EOFir ingPr ess . Text Then 
EarlyEnd = i 
Exit For 
L End If 

Next i 

InPoints = BarlyEnd 
-If Sgn (Goal - ThisValue) <> Sgn (Direction) Then 
Step = Step / 2 

Direction = Sgn (Goal - ThisValue) 
End If 

plf Sgn (Goal - ThisValue) > 0 Then 

a = a + Step 
-Else 

a = a - Step 
Lend if 

-Loop until .Abs(Goal - ThisValue) <= Accuracy 

StatusLabel. Visible = False 
StatusText .Visible = False 
ProgressBarl .Visible = False 
DataCut Frame. Visible = True 
Refresh 

End Sub ^ 

Private Sub ThrustCoef Pactor_MouseDowrL (Button As Integer, Shift As Integer, X 
As Single, Y As Single) 

cfFactorPopup. visible = True 

End Sub 



Private Sub ThrustCoef Factor_MouseUp (Button As Integer, Shift As Integer, X 
As Single, Y As Single} 

CfFactorPopup .Visible = False 

End Sub 



FIG. 3B 



Private Sub SetBurnRate () 

'This subroutine varies the reference bum rate at the input temperature 
'until the web burned at end of firing calculated from the input pressure 
'data matches the theoretical (or user entered) web. This subroutine 
'builds the web burned array, which becomes part of the burn back file 

Dim ThisValue As Double, Goal As Double 

Dim step As Double, Direction As Integer, Accuracy As Double 




Dim a As Double, b As Double, c As Double, d As Double 

Dim rref As Double, brate As Double, trefin As Single, trefout As Single 
Dim nexp As Single, PrevBr As Double, Trial As Integer 

DataCutFrame. Visible = False 
StatusLabel. Visible = True 

StatusText .Caption *= "Calculating Reference Burn Rate..." 
StatusText .Visible = True 
ProgreasBarl .Value = 0 
ProgressBarl. visible = True 

Goal = Val (inWebDepth. Text) 
Accuracy = 0.001 

Direction = 0 'first time thru step will be halved 

Step = 0.1 'before it is applied 

rref = Val (InBurnRate .Text) 

trefin = Val (InTemp .Text) 

nexp = Val (InN. Text) 

PrevBr = Val (InBurnRate .Text) 

Trial = 0 

-Do 

Trial = Trial +1 

Thisvalue = o 

For i ■ l To inPoints 

plf 1 / int (InPoints. / 10) = lnt(i / lnt(inPoints / 10) ) Then 

StatusText .Caption = "Calculating Reference Burn Rate... Trial #» & Trial 
ProgressBarl. value = i / InPoints * loo 
Refresh 
End If 

j brate = rref * (Pressln(i) / PREF) A nexp 'eg 4 

Thisvalue = Thisvalue + (PrevBr + brate) / 2 * (Timeln(i) - TimelnU - l)) 
'eg 5 

WBumln(i) = ThisValue 
PrevBr - brate 
Next i 

plf Sgn(Goal - Thisvalue) <> Sgn (Direction) Then 
step = step / 2 

Direction = Sgn (Goal - Thisvalue) 
L~End If 

If sgn (Goal - Thisvalue) > 0 Then 

rref • rref + step 
-Else 

rref = rref - Step 
•-End If 

Loop Until Abs (Goal - ThisValue) <= Accuracy 

StatusLabel .Visible « False 
StatusText .Visible » False 
ProgressBarl. Visible = False 
DataCutFrame. Visible = True 
Refresh 

InBurnRate .Text = rref 

OutBumRate.Text = Val (InBurnRate .Text) * Exp { (Val (SigmaP. Text) * (Val (OutTemp. 
Text) - Val (InTemp. Text) ) ) ) 

End Sub 



FIG. 4B 



Privat Sub SstCStarO 

This function varies the value of C* until the propellant weight burned 
'at the end of firing matches the given propellant weight. This routine 
'builds the burn surface and thrust coefficient arrays that become part 
'of the burn back file. 

Dim ThisValue As Double, Goal As Double 

Dim Step As Double, Direction As Integer, Accuracy As Double 
Dim a As Double, b As Double, c As Double, d As Double 

Dim rref As Double, brate As Double, trefin As Single, trefout As Single 
Dim nexp As Single, mpdot As Double, PrevMpdot As Double, CS As Double 
Dim rho As Double, Trial As Integer 
Dim bsmax As Double 



Da taCut Frame .Visible = False 
StatusLabel .Visible = True 

StatusText .Caption = "Matching Propellant Weight..." 
StatusText .Visible ■ True 
ProgressBarl .Value = 0 
ProgressBarl .Visible = True 

Goal = Val (Pwt.Text) 
Accuracy = 0 .01 

Direction - 0 'first time thru step will be halved 

Step = 1000 'before it is applied 

rho = Val (InDensity. Text) 

rref = Val (InBurnRate. Text) 

trefin = Val (InTemp .Text) 

nexp = Val (InN. Text) 

cstar.Text = eoooo / 12 

CS = Val (CStar.Text) * 12 

Trial =0 

-DO 

Trial = Trial + l 
Thisvalue - 0 
bsmax =0 

For i = 1 To InPoints 

r lf i / mtdnPoints / 10) = lnt(i / intdnPoints / 10}) Then 

StatusText .Caption - "Matching Propellant Weight... » & Chr$(l3) & "Trial 
#" & Trial " * 

ProgressBarl .Value = i / InPoints * 100 
Refresh 
L End If 

mpdot ■ Pressln(i) * Atln(i) * G / CS 'eg 6 
brate = rref * (Pressln(i) / PREF) a nexp 'eg 4 
BSurfln(i) = mpdot / (brate * rho) » e q 4 

If bsmax < BSurfln(i) Then 
bsmax = BSurfln(i) 
L End If 

Thisvalue = ThisValue + (mpdot + PrevMpdot) / 2 * (Timeln(i) - TimelnU - 1) ) n 
'eq 8) v 
PrevMpdot = mpdot 

cfln(i) = Thru3tln(i) / (Pressln(i) * Atln(i)) *eq 9 
L-Next i 

If Sgn(Goal - ThisValue) <> Sgn (Direction) Then 
Step = step / 2 

Direction = Sgn (Goal - ThisValue) 
^End if 

plf Sgn (Goal - ThisValue) > 0 Then FIG 5 A 



cs = cs - step 
—Else 

cs . cs + step FIG. 5B 

End If 

-Loop Until Abs (Goal - ThisValue) <= Accuracy 

StatusLabel .Visible = False 
StatusText .visible = False 
ProgressBarl .Visible = False 
DataCutFrame. Visible = True 
Refresh 

CStar . Text = Format (CS / 12 , 11 o 11 ) 
MaxBS = bsraax * 1.2 



find Sub 



Private Sub SetErosion () 

'This subroutine varies the coefficient a in the equation E=aP*b 

'until the final throat calculated matches the given throat. 

'E is the incremental radial throat erosion and P is the incremental 

'chamber pressure. The coefficent b is fixed at .924. This routine build the 

throat area area that becomes 

'part of the burn back file 

Dim ThisValue As Double, Goal As Double 

Dim Step As Double, Direction As Integer, Accuracy As Double 

Dim a As Double, b As Double, c As Double, d As Double, Trial As Integer 

DataCutFrame. Visible = False 
StatusLabel .Visible = True 

StatusText .Caption = "Calculating Erosion Coefficients • » • " 

StatusText .Visible «= True 

ProgressBarl. Value - 0 

ProgressBarl .Visible = True 

Refresh 

'Calculate erosion coefficients 

a = 0.00001477 
b = 0.324 

Goal = Val (AtF.Text) 
Accuracy ■ o.oooooooi 

Direction = 0 'first time thru step will be halved 

Step = o.ooooi 'before it is applied 

Trial = 0 



-Do 

Trial = Trial + l 
ThisValue = Val (AtO .Text ) 
EarlyEnd = InPoints 
For i o i to InPoints 
r If i/ Int (InPoints / 10) = Int(i / Int (InPoints / 10)) Then 

StatusText. Caption = "Calculating Erosion Coefficients... Trial # n & Trial 
ProgressBarl .Value = i / InPoints * 100 
Refresh 
End If 

ThisValue = ThisValue + (a * Pressin(i) A b) *. (Timeln(i) - TimelnU - 1)) * 
2 'eg (2 J 

Atln(i) = ThisValue A 2 / 4 * PI 'eq(3) 



Private Function Oener ateOutput ( ) 

'This function uses the Jbura surface profile, cf profile and throat area 
'pxvfile generated from the input data to generate a prediction at the 
'output temperature. The model usee an iterative proceee that varies a 
'bum surface modifier until the calculated propellant weight burned matches 
'the actual propellant weight. 

Dim Thisvalue As Double, Goal As Double 

Dim Step As Double, Direction As Integer, Accuracy As Double 
Dim a As Double, b As Double, c As Double, d As Double 

Dim rref As Double, brate As Double, trefin As single, trefout As single 
Dim nexp As Single, mpdot As Double, PrevMpdot As Double, cs As Double 
• Dim rho As Double, Trial As Integer 

Dim BSFactor As Double, BSModifier As Double, PercentWeb As Single 

DataCutFrame. visible - False 
StatusLabel. visible = True 

S tatusT ext. Caption = "Generating Converted File... 1 ' 
StatusText .Visible = True 
ProgressBarl .Value = 0 
ProgressBarl .Visible = True 

Goal = Val (Pw tout. Text) 
Accuracy = 0.05 

Direction = 0 'firet time thru step will be halved 

Step =0.02 'before it is applied 

rho o val (OutDensity. Text) 

rref = Val (OutBurnRate .Text) 

trefin = Val (InTemp. Text) 

trefout = Val (OutTemp. Text) 

PressOut(Q) =14.7 

nexp = val (OutN. Text) 

CS = val (est ar. Text) * 12 

BSModifier = 1 




Trial b o 



-Do 

Trial = Trial + 1 
Thisvalue = o 
For i = 1 To InPoints 
plf i / int (InPoints / 10) = Int(i / Int (InPoints / 10)) Then 

StatusText .Caption = "Generating Converted File... " & Chr$(13) & "Trial #" 

& Trial 

ProgressBarl .Value = i / InPoints * 100 
Refresh 
L-End if 

WBurnOut(i) = WBurnln(i) * OutWebDepth . Text / InWebDepth . Text 'eg 10 
Percentweb = WBurnOut { i ) / OutWebDepth . Text 
BSurfOut(i) a BSurfln(i) * BSModifier 'eg 11 

brate » rref * (PressOut (i - 1) / pref) a nexp » e g 12 

mpdot = brate * (BSurfOut (i) + BSurfOut (i -1)) / 2 * rho 'eg 13 
plf AtFOut.Text <> AtO.Text Then 

AtOut(i) = (AtO.Text + (Sqr (4 * Atln(i) / PI) - AtO.Text) / (At P. Text - AtO 
.Text) * (AtFOut.Text - AtO.Text)) A 2 / 4 * PI 'eg 14 
-Else 

AtOut(i) = Atln(i) 
End If 

' Values from form (AtO.Text, etc) are actually diameters, not areas) 
I PressOut (i) = (mpdot * CS) / (AtOut(i) * G) 'eg 15 
If PressOut (i) < 14.7 Then 

PressOut (i) =14.7 
End If 

CfOut(i) = cfln(i) * (l + (ThrustCoef Factor .Text / 100) * (trefout - trefin) 
/ 100) 'eg 16 

ThrustOut(i) = Cfout(i) * PressOut (i) * AtOut(i) 'eg 17 

TimeOut (i) = TimeOut (i - 1) + (WBurnOut (i) - WBurnOut (i - 1) ) / brate 'eg 
18 

Thisvalue = ThisValue + (mpdot + PrevMpdot) / 2 * (TimeOut (i) - TimeOut (i - 1 
)) 'egW. 
PrevMpdot - mpdot 

BBPGraph.Llne (WBurnOut (i - l) , MaxBS * 1.3 - BSurfOut (i - l ))- (WBurnOut ( i ) , 
MaxBS * 1.3 - BSurfOut (i) ) , BLUE 

ThrustGraph.Line (TimeOut (i - 1), MaxThrust * 1.3 - ThrustOut(i - 1) ) - 
(Timeout (i) f MaxThrust * 1.3 - ThrxistOut ( i ) ) , BLUE 

Next i 

-If Sgn(Goal - ThisValue) <> Sgn (Direction) Then 
Step «* Step / 2 

Direction = Sgn (Goal - Thisvalue) 
End If 

If Sgn (Goal - ThisValue) > o Then 
BSModifier = BSModifier + Step 
-Else 

BSModifier = BSModifier - Step 
End If 

-Loop until Abs (Goal - Thisvalue) <= Accuracy 

StatuaLabel .Visible = False 
StatusText .Visible = False 
ProgressBarl. Visible = False 
DataCutFrame. Visible - True 
Refresh 

CS tar. Text = Format (CS / 12 , "0") 
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